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13  ABSTRACT 


A  theory  is  derived  for  the  flexural  deformation  of  laminated  beams  subjected 
to  initial  stresses.  Each  layer  is  treated  as  a  Timoshenko  beam  and  "smoothed" 
expressions  for  the  kinetic  and  strain  energies  and  the  work  do.ie  by  external  forces 
are  derived.  A  system  of  three  coupled  partial  differential  equations,  including 
the  effect  of  bending,  extension,  rotation,  and  initial  axial  stress  are  obtained 
from  Hamilton's  principle.  Specific  buckling  and  free  vibration  problems  are  solved 
exactly  for  hinged-hinged  and  clamped-clamped  beams.  The  numerical  results  reveal 
values  of  buckling  coefficients  and  natural  frequencies  that  are  in  agreement  with 
the  corresponding  results  obtained  from  the  effective*" modulus  theory  fo*  relatively 
long  beams,  but  that  are  considerably  lower  than  the  effective  modulus  values  for 
relatively  short  beams.  In  the  high  frequency  range,  with  the  effect  becoming  more 
pronounced  in  the  higher  modes,  the  present  theory  predicts  much  smaller  values  for 
the  flexural  and  thickness-shear  natural  frequencies  than  does  the  effective  modulus 
theory. 
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ABSTRACT 

A  theory  is  derived  for  the  flexural 
deformation  of  laminated  beams  subjected 
to  initial  stresses.  Each  layer  is  treated 
as  a  Timoshenko  beam  and  "smoothed**  expres¬ 
sions  for  the  kinetic  and  strain  energies 
and  the  work  done  by  external  forces  are 
derived.  A  system  of  three  coupled  partial 
differential  equations,  including  the  effect 
of  bending,  extension,  rotation,  aiid  initial 
axial  stress  are  obtained  from  Hamilton's 
principle.  Specific  buckling  and  free 
vibration  problems  are  solved  exactly  for 
hinged-hinged  and  clamped-clamped  beams. 

The  numerical  results  reveal  values  of  buck¬ 
ling  coefficients  and  natural  frequencies 
that  are  in  agreement  with  the  corresponding 
results  obtained  from  the  effective  modulus 
theory  for  relatively  long  beams,  but  that 
are  considerably  lower  than  the  effective 
modulus  values  for  relatively  short  beams. 

In  the  high  frequency  range,  with  the  effect 
becoming  more  pronounced  in  the  higher  modes, 
the  present  theory  predicts  much  smaller 
values  for  the  flexural  and  thickness-shear 
natural  frequencies  than  does  the  effective 
modulus  theory. 
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INTRODUCTION 


Recently  C.-T.  Sun  [1]  developed  a  natheiaatical  theory  for  a 
continuum  model  with  microstruccure  for  laminated  beams.  These 
laminated  beams  were  assumed  to  cor^vist  of  many  parallel  alternat¬ 
ing  layers  of  two  homogeneous,  isotropic  elastic  materials.  Each 
constituent  layer  was  treated  as  a  Timoshenko  beam,  and  expressions 
for  the  kinetic  and  strain  energies,  as  well  as  the  work  done  by 
external  forces,  were  derived.  The  smoothing  technique,  which  had 
previously  been  used  so  effectively  in  deriving  a  general  field 
theory  for  laminated  composites  [2-5],  was  applied,  and  the  equa¬ 
tions  of  motion  and  associated  boundary  conditions  were  derived  from 
Hamilton’s  principle.  To  test  the  accuracy  of  the  resulting  theory, 
the  dispersion  curve  for  flexural  waves  in  an  infinite  laminated 
beam  were  determined  from  (i)  the  exact  theory,  (ii)  the  new  theory, 
and  (iii)  the  effective  modulus  theory  for  a  composite  beam  consisting 
of  five  stiff  layers  and  four  soft  layers.  It  was  found  that  the 
curves  obtained  from  the  exact  theory  and  the  microstructure  theory 
were  in  excellent  agreement  over  the  entire  range  of  wave-number 
considered.  By  contrast,  the  effective  modulus  theory  exhibited 
agreement  only  in  the  very  low  frequency  range. 

The  effective  modul\is  theory  treats  the  laminated  composito  as 
a  homogeneous  but  transversely  isotropic  medium.  Based  on  this  point 
of  view,  Brunei le  [6*5]  has  investigated  the  stability  and  vibration 
characteristics  of  transversely  isotropic  Timoshenko  beams  imcler 
initial  stress.  He  showed  that  the  values  of  buckling  loads  and 
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natural  frequencies  of  vibration  are  strongly  dependent  upon  the  nature 
of  the  boundary  restraint  and  the  relative  dejjree  of  anisotropy,  i.e,, 
the  ratio  of  the  longitudinal  Young's  modulus  to  the  transverse  shear 
modulus.  Decreases  in  the  values  of  the  buckling  coefficients  and 
natural  frequencies  of  vibration  relative  to  the  classical  values, 
computed  on  the  basis  of  the  Euler-Bemoulli  theory  for  beams,  were 
observed  to  be  particulaxly  pronounced  in  the  case  of  the  clamped- 
clamped  beam. 

The  objective  of  the  present  investigation  i'j  to  extend  Sun's 
microstructure  theory  for  laminated  beams  [1]  so  as  to  include  the 
effect  of  initial  stresses  and  to  solve  a  pair  of  stability  find  free 
vibration  j’voblems  within  the  frsmevv-'jrk  of  the  resulting  theory. 

The  numerical  results  for  buckii  fy  coefficients  and  frequency  para¬ 
meters  are  compared  with  Br?»r;;iie's  results  [6-9]  derived  from 
Timoshenko  beam  theory  for  transversely  isotropic  materials. 


STATEMENT  OF  THE  GENERAL  PROBLEM 


We  consider  an  elastic  solid  of  volume  V  bounded  by  a  finite 
surface  S.  Displacements  are  prescribed  on  the  portion  S^  of  the 
surface  S  and  tractions  are  applied  on  the  remaining  portion  S^, 
Initially  the  body  is  assumed  to  be  at  rest  and  is  subjected  to  a 
state  of  initi.'..'  stress  i,  j  a  1,2,3,  which  arises  from  con¬ 
servative  forces  applied  on  the  surface  S^  of  V,  We  assume  that 
the  perturbed  aiid  linearized  field  equations  [10,11]  are 


'ij,j  *  '^“jk  *  ''i  “  ““i- 
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t'ij  *  “jk  “i,k)”j  ■  ™ 

prescribed,  on 


‘^ij  “  %»■* 


T 


(2) 

(3) 

(4) 


where  denotes  the  perturbed  state  of  stress,  is  the  displace¬ 
ment  v«£Ctor  measured  from  the  undisturbed  state.  is  the  body  force 
field  per  unit  volume,  p  is  the  mass  density,  and  n^  is  the  unit 
exterior  normal  to  the  surface  S,  The  quantity  6  is  a  parameter 
associated  with  the  magnitude  of  the  externally  applied  surface 
tractions,  and  the  pj  are  the  components  of  perturbations  of  the 
applied  surface  tractions  on  In  (4),  the  general  stress-strain 
relationship  is  given.  The  familiar  summation  convention  for 
Cartesian  tensors  is  employed  in  (l)-(4),  commas  denote  spatial 
derivatives,  and  dots  time  derivatives. 


THE  VARIATIONAL  FORMULATION 

In  order  to  derive  a  form  of  Hamilton's  principle  that  will 
generate  (l)-(3),  we  multiply  (1)  by  Su^  and  form  the  integral  of 
the  volume  V  of  the  result: 


/  T..  .6UkdV  +  6/  .)  .6u.dV  +  /  F.du.dV  ■  /pu.i^u.dV. 

Y  lj,J  y  tJ-*- 

(5) 

Application  of  the  divergence  theorem  to  (5)  yields 

I  (^ij  ♦  Bajj^Ui^^)nj6UidS  -  /6(W  *  j  Bo * 
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RjBMMl»irSuqa«n 


+  /  F.Su.dV  =  /  pu.6u.dv, 

J  1  1  ^  1  ^ 


where  we  have  exploited  the  symmetry  properties  of  and 
to  write 

6W  =  (3IV/3e..)6e..  =  T..6e..  =  t..6u.^. 

^jk“i,k'^”i,j  "  7  ^^°jk“i,k“i,j^* 

with  W  denotinn  the  strain  energy  density.  Furthermore,  if  we 
assume  that  the  body  force  vector  can  be  represented  as 


(6) 


=  ^i  ^  ^i j  ~ 


then 


F  6u.  =  '5(f  u.  +  i-  a.  .u.u.), 

i  1  i  1  2  13  1  3 


Hence,  in  view  of  (2)  and  (3),  (6)  can  be  expre-:sed  as 

(5  /  [W  +  ~  Bo.,  u.  ,  u.  .  -  f.u.  -  ~  ai^UiU.jdV  + 

i,  ^  2  ]k  i,k  i,j  1  1  2  13  1  y 


+  /  pu.6u, dV  =  6/  Bp.u.dS, 

y  1  1  S  ^  ^ 


(7) 


Forming  the  time  integral  of  (7)  and  integrating  by  parts  in  the 
term  involving  u^  6uj,  wc  obtain 

<5  /  Xdt  =  0, 

t„ 


(8) 
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where 


;£  ■  /  [  i  PVi  -  T  -  *l"i  -  T  * 

♦  /  Sp.u.dS. 

S,  P) 

Eq.  (8)  represents  Hamilton's  principle  associated  with  the 
problem  (1)-C3),  In  conjunction  with  (9),  (8)  can  be  used 
effectively  for  the  purpose  of  deriving  plate  and  beam  theories 
when  initial  stresses  are  present. 


SOME  PRELIMINARY  CONSIDERATIONS  FROM  BEAM  THEORY 

To  derive  an  elementary  theory  for  extensional  deformations 
and  a  Timoshenko  theory  for  flexural  deformations  in  beams,  we 
approximate  the  displacement  field  as  follows: 

u^(x,t)  =  u(xpt)  -  X2<KXj,t),  U2(x,t)  »  w(xi,t) 

u,  (x,t)  *  0, 

(10) 

where  the  Xj  -  axis  coincides  with  the  axis  of  the  beam  which 
passes  through  the  centroid  of  every  cross-section  and  the  X2  - 
axis  denotes  the  transverse  direction.  Here  u  is  the  longitudinal 
deflection  of  the  axis  of  the  beam,  w  the  transverse  deflection, 
and  <Ji  the  rotation  of  a  section  in  the  XjX2  -  plane.  The  strains 
computed  from  (10)  are  easily  shown  to  be 

hi  =  '".l  -  V,1  h2  ”  J  (”,1  -  h3  “  hs  “  ^22  "  "33 

(11) 
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In  the  usual  manner,  we  shall  ignore  the  relations  »  0  in 

favor  of  setting,  for  an  isotropic  material, 

T22  =  (2u  ♦  X)e22  ♦  ^(Ejj  ♦  £33)  ■  0, 

"^33  ~  *  ^^^33  ■*■  ^^^11  *  ^22^  “ 

Thus,  it  follows  that 

^22  *  *^33  “  "  *  W),  »  UEjj/Cu  +  A). 

(12) 

The  expression  for  the  strain  energy  is 


ueijEij 


=  liCe^j  +  C22  ♦  elj)  ♦  2u(Cj2  +  Cjj  ♦  e^j)  +  i-Ae^^c^^. 

(13) 


Substitution  of  (11)  and  (12)  into  (13)  leads  to 


(U) 


since  E  «  y (2u+3A)/(u+l) ,  E  being  Young's  modulus.  In  the  development 

of  the  Timoshenko  beam  theory,  it  is  customary  to  introduce  a  shear 

2  2 

correction  factor  <  into  (14),  Hence,  we  replace  £^2  by  <£22* 


(14)  becomes 


W  =  i  Ee^j  +  2yKe22, 


or,  by  virtue  of  (11) 

W  =  IfCu'i  -  .  x2«2j)  .  ixuCw^l  - 


(15) 
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An  integration  of  (15)  over  the  cross-sectional  area  A  of  the  beam 
yields 

U  *  /  WdA  «  i  EAu^,  ♦  I*  ^  kuA(w  ,  - 

^  2  Z  ,1  2  ,1 

(16) 

where  I  »  /  x^dA  and  /  x.dA  »  0  since  the  x.  -  axis  passes  through 
the  centroid  of  the  cross-section. 

In  a  similar  fashion,  for  the  kinetic  energy  density,  we  have, 
in  view  of  (10), 

0  ‘  J  ^  P(“^  “  2*2“^  *  ^2^^^  *  jPOf^ 

and  thus 


T  »  /  i.  ou.u.dA  »  i  p(Au^  +  1$^  *  Aw‘), 

A  2  1  X  2 


Furthermore,  we  can  also  show  that 


(17) 


'’jk“i,k“i,j  ■  '’at'*!,!  *  “2,1  *  “3,1’  *  ‘’22t'‘l,2  *  "1.2  *  “I,2>  * 


*  '’33(“l3  ♦  “2,3  *  'i.3>  *  2°12!"i,2“i,1  * 


*  “2,2“2,1  *  “3,2“3,!'  *  ^‘'l3t“l,3"l,l  *  “2,3“2,1  * 


*  “3,3“3,1>  ♦  2'’23f“l,3"l,2  *  “2,3“2,2  *  "3,3“3,2’- 

If  we  assume  that  all  »  0  except  o^j,  then,  in  view  of  (10),  we 
have 


10 


and 


If  »  0,  then,  using  (16)- (18),  we  can  reduce 


(9)  to 


=  /  (T  -  U  -  U*)dx-. 


ENERGIES  FOR  A  LAMINATED  BEAM 


In  this  section  we  intend  to  extend  the  theory  developed  by 
Sun  fl]  for  laminated  composite  beams  to  include  the  effect  of  the 
initial  stress  From  this  theory  we  should  be  able  to  compute 

buckling  loads  for  such  beams.  We  consider  a  composite  beam  con¬ 
sisting  of  a  large  number  of  parallel  alternating  layers  of  two 
homogeneous,  isotropic  elastic  materials  (See  Figure  1).  We  assume 
that  the  cross  section  of  each  layer  is  rectangular  and  of  width  b. 

The  depth  of  the  cross-section  shall  be  denoted  by  h.  We  shall 
designate  the  properties  of  the  stiff  layer  by  the  subscript  1  and 
of  the  soft  layer  by  the  subscript  2.  The  thicknesses  of  the  layers 
are  dj  and  ^2*  Young's  moduli,  shear  moduli,  and  the  mass 

densities  for  the  two  types  of  constituents  are  denoted  by  Ej,  Uj, 
and  E2,  U2,  P2»  respectively.  In  addition,  we  shall  use  a  superscript 
k  to  refer  to  the  properties  of  the  k-th  stiff  and  soft  layers. 
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According  to  (16),  the  strain  energies  per  unit  length  in  the 
k-th  stiff  and  soft  layers  due  to  bending,  shear  deformation,  and 
extension  are 


“5  -  T  '^1*1  ("5.1)^*  2  T  "I"!*!'"?,!  - 


k,2 

» 


”2  “  }  l^^""  2'''2''2'^2f”2  1  ’  * 


respectively,  where,  as  is  evident  from  Figure  1, 


Aa  =*  ^  =  K/12.  «  =  1,2. 


By  virtue  of  (18),  the  potential  energies  per  unit  length  due 
to  the  initial  stress  are 


Similarly,  the  kinetic  energies  per  unit  length  for  the  k-th 
stiff  and  soft  layers  are 

=  i  OjIAjCu^)^  *  ‘ 


"M-%'"2'“2’'  *  ■2t*2>'  *  ''2'V)'1- 


resnectively,  in  view  of  (17). 
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Sun  [1]  introduced  an  additional  approximation  on  the  extensional 
components  namely. 


Uq  =  -  a  »  1,2, 


(27) 


where  '^'Cx^jt)  denotes  the  gross  rotation  of  the  composite  beam  and 

V 

y^  the  coordinates  of  the  centroidal  axes  of  the  k-th  stiff  and  soft 
layers.  Insertion  of  (27)  into  (20) -(26)  leads  to 

*  f  *  T 

C28) 

C  ■  ?°n  A0'a)%^  ‘ 

(29) 


(30) 


Onder  the  assumption  that  there  are  n  pairs  of  stiff  and  soft 
layers  in  the  composite  beam,  the  total  strain,  kinetic,  and  potential 
energies  per  unit  length  can  be  expressed  as 


(U,  T,  U*)  -  I  f  (uj,  tJ,  U^;*'). 

k*l  0*1 

(31) 

Still  following  Sun  [1],  we  introduce  a  smoothing  operation  through 
which  the  summation  over  k  in  (31)  may  be  approximated  by  the 


following  weighted  integrations  over  y; 


(U,  T,  U*)  »  / 


■h/2  1 


-h/2  dj+d2 


1. 


t»l 


(32) 
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Ic 

where  and  y^  are  to  be  re^ laced  by  y  and  the  superscript  k  is 
deleted^  Eq,  (32)  is  also  valid  for  the  crss  in  which  the  number 
of  stiff  layers  is  not  equal  to  the  number  of  soft  layers. 

Inserting  C28-(30)  into  (32),  we  obtain 
2 

^  2 


U*=i 


where  C  =»  h/(dj+d2). 

The  number  of  dependent  variables  appearing  in  (33)- (35)  can  be 
reduced  because  the  continuity  conditions  at  the  interface  between 
stiff  and  soft  layers  introduces  two  relationships  involving  w^, 
and  ij>.  In  terms  of  the  disci'st©  variables,  the  continuity  of 
displacements  at  the  intorfacs  of  the  k-th  pair  of  layers  leads  to 


“■  -  "2  ““  "i  *  ? ’’ih  =  "2  -  f  V2 


'  r  k 

^  T  I  d 


But  in  view  of  (27),  (36)  becomes 


Cyk  .  .  i  ^ 

-I  ^  ^  OMl  ° 

and  since  y ,  ®  Xo  ♦  T  result  reduces  to 

>}»  a  n<{ij  +  (l-n)*}*^* 

where 

n  «  dj/ (d^+d^) . 


(37) 


(38) 


Applying  the  smoothing  operation  to  (37),  we  have 


li*  ■  nijij  +  (l-n)<l'2» 


Solving  for  we  find 


41^  a  (ti(-n'(')/(l-n), 


(39) 


where  we  now  set  =  <l>.  Therefore,  substitution  of  (39)  and  Wj 


Wj  =  w  into  (33)- (35)  yields 


E  I 

*  I  ^  i'P  1^'*, 


2  '  (i_n)2  2  '^WZ  }  ",1 


illllliill 

1-n  J 


(40) 


£  Oi,I.. 


*1  2l  21  2^  ”11*2  2 

'J  *■  I"*"  1  ■*’  . “""T  P  * 

2  11  b  ,1  Z  11  2  11  1  ,1  2(1-1))^  * 


(41) 
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(44) 


.  "  . . .  .  . . ,  , ^  ■^.  ,  -  .  ■  ■  -  -  -  -  . . - ■■.-■■  

^3  “  *  icbdjCiii  "  1*2^ 

^4  “  *^1^1  *  ^2^1  “  *  ^2^2^' 

ag  »  E2l2/Cl-n)^  +  +  (l-n)E2.y5 

-  (b/12)  [E2d2(dj^  ♦  d2)2  +  h^CEjdi  +  E2d2)l. 

=  a2/(l-n)  "  (b/d2)tcy2  (di  +  d2)‘‘it 
ay  a  nE2l2/(l*‘’l)^  *  C>/12)E2djd2Cij  ♦  d2)» 
ag  *  nag  »  Cb/d2)KP2‘^l  ^*^1  '^2)» 

ig  *  p2l2/Cl-n)^  *  np^ys  ♦  (l“n)p2V5 

«  Cb/12)  [P2d2  (dj  +  d2)^  ♦  h^(Pjdj  ♦  P2d2)l » 

*10  “  nP2l2/Cl”’^)^  *  Cb/12)p2did2  (dj  +  dy), 
ail  "  h^l*  n^E2l2/(l"n)^  -  (b/12)df  (Eidi  +  E2d2). 

®12  “  ■  Cb/d2)tcdi(Wid2  +  y2^l)» 

*13  “  *  Cb/12)di  (Pidi  +  P2‘J2^» 

ai4  -  ■*•  X^/Cl-n)^  =  Cb/12)(di  +  d2)[h2  +  d2(di  +  d2)'j , 

^IS  “  h  *  “  Cb/12)d^(di  +  y, 

(45) 

where  a  ^2  =  k  »  0,822  for  a  rectangular  cross-section. 
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Insertion  of  (44)  into  Hamilton's  principle  (8)  leads  to  the 


following  system  of  displacement  equations  of  motion  in  0<Xj<t: 
(a^  -  a2<J'^i  -  »  a^w. 


+  ag^.  »  agiji  -  a^Q^., 


®3”,i "  ®ii/^2^'^,n  *  ^®ii  ■*■  ®15°ll^‘*’,n  “ 


-  aj2<f>  **  a^jv  -  ajQi|<, 


with  the  boundary  conditions 


(i)  either  (a^  +  ajjA/5)w^l  -  a2'J»  -  a.^<^  =  0  or  5w  «  0, 

(ii)  either  (Ag  +  aj^o^j)!!)^!  -  a^Cl  +  on/E2)ii)^  >  *  0  or  6iJ)=  0 

(iii)  either  (ajj+  i  ”  ^7^^  ■*■  015/^2^’*'  1  ~  ^  or  6(ji*  0 

(49) 

on  Xj  a  0,5.,  If  we  set  Ojj  »  0  in  (46)-(49},  then  Inese  equations 
become  identical  to  Sun's  equations  (45)-(50)  in  [1]  in  the  absence 
of  surface  and  end  loads. 

Let  us  now  suppose  that  a  compressive  load  P  is  applied  along  the 
axis  of  the  beam,  so  that  Oj^jA  =  -  P,  Then  (4L)-(48)  become 

(bj-Xb2)w"  (x,T)-lir' (x,T)-.b3?' (x,t)  »  bj3W(x,T), 


9 


(d^/£)^(Eidi>E^d^)  6Cc6j)2(1^ 


12KUzCdj+dz) 

$ 

12ic(l+d) 

Ez(di/£)2CiTh/£)2 

144^^2 

144ic 

diCuidz+Uzdi) 

d(Y+d) 

1+d 

1  +  0d 

icp^Cdj+dz) 

K(l-*-d) 

p2d2(tli+dz)^+h^(pldi-*-P2d2)  [B^  (l+d)^+l+0d] 

12K£^P2Cdj+d2)  12ic(l+d) 

djdj  ?2Sj8^ 

12k£^  12k 

$ 

^1 1^1^^2^2^ 

6j;^(i+ed) 

12K£^Pz(dj+d2) 

12k  (l+d) 

Ej/Ez,  y 


Uj/U2.  ^ 


dj/dz. 


C  =  h/£, 


pjp^ 


fi 

a 


d  /h, 

a 


o  =»  1,2, 
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SOLUTION  OF  TliE  FIELD  SQUA.T10NS 
The  Hinged-Hinged  Beam 

We  now  consider  the  motion  of  a  beam  that  is  hinged  at  the  ends 
X  =  0,1,  According  to  (49),  the  boundary  conditions  are 

w  -  (b,^-Xb5)?•  -  b^(l-Xbg)?'  »  b7Cl-Xbg)?  -  Cbio-^bjj)7'  =  0 


(5. -5) 


at  X  a  0,1.  It  is  a  simple  matter  to  verify  that  the  functions 


w(x,t)  s>  A2  sin  (mix)  cos  wt,  Hx,t)  *  A2  cos  (mix)  cos  wx, 


<|i(x,t)  =  Aj  cos  (mix)  cos  wr,  n  »  1,2,3,,,., 


(54) 


satisfy  the  boundary  conditions  (53)  identically.  Substitution  of 
(54)  into  the  equations  of  motion  (50)- (52)  yields  the  homogeneous 
algebraic  system 


Oij  Aj  •  0, 


(55) 


where 


°11  *’13“^  ■  “22  “  ’^14“^  ^  ^6  ”  Cn7i)^(b4-Xb5), 


2  7 

“53  "  *^16“  '•  ^12  “  ^lO'^'^lP'  “19  ”  “ 


12  21 


IITT 


2  2 

“23  =  “32  °  b7(l-Xbg)  -  bj^u  ,  =  niib^. 


(56) 


22 


If  (55)  is  to  have  a  nontrivial  solution,  we  must  require  that 


det  (oj^.)  a  0. 

(57) 

Tlie  values  of  the  natural  frequencies  and  of  the  critical  load 
parameter  (o)=0)  will  be  computed  from  (57). 


The  Clamped-C lamped  Beam 

For  combinations  of  boundary  conditions  other  than  those  of  the 
hinged-hinged  c'_-e  described  above,  simple  expressions  of  the  form 
(54)  for  the  solution  of  (50)- (52)  cannot  be  found.  Nonetheless,  we 
can  obtain  the  values  of  the  natural  frequencies  and  critical  loads 
for  these  more  difficult  cases  by  following  a  procedure  described  by 
Jones  [12].  If  we  insert 

w(X,t)  a  w(x)  cos  OJT,  T(X,T)  =  i;»(x)  COS  UT, 

^(x,T)  =  I^x)  cos  UT, 

into  (50)- (52),  we  find  the  following  set  of  ordinary  differential 
equations; 

C  W»*(x)  +  C,w(x)  -  l(»'(x)  -  C^ij>'(x)  a  0, 
i  ^  4 

W'(X)  +  C5lJj*'(x)  +  Cgl|;(x)  -  Cy<j)”(x)  +  Cg<|i(x)  a  0, 

C4W'(X)  -  Cyi|<'’(x)  +  Cg'Hx)  +  C9i|>"(x)  <•  Cj()Mx)  =  0, 

(58) 
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1 


where 


Cl  «  C2  •  C4  -  bj,  Cg  »  b4-Ab5, 

Cg  »  bi4aj2»bg,  «  bTd-Xbg),  Cg  »  bg-bi3c 


2 

;o)  , 


Cg  »  bio-Xbii,  Cjg  »  bi^a>2  .  bi2. 


(59) 


For  the  case  of  buckling,  we  set  u-0,  so  that  (58)  reduces  to 


Cj^w“(x)  -  ij)'(x)  -  »  0, 


w'(x)  +  Cgij/'Ux)  ♦  Cg'i'(x)  -  C2'l'**(x)  +  Cg(ji(x)  ■  0, 


C4w'(x)  -  C2<<»"(x)  ♦  Cgil'(x)  Cg'}'**(x)  ♦  Ciq()>(x)  ■  0, 


(60) 


where  now 


S  “  0»  ‘'8  ’  ^9»  ‘'lO  • 


^12’ 


Since  the  equations  in  (58)  and  (60)  possess  constant 
;oofficients,  we  seek  a  solution  in  t!'o  form 


w(x)  ■*  Ae^^,  "{((x)  *  4-(x)  =  Ce'^, 


Substitution  of  (61)  into  (58)  yields  the  following  homogeneous 
system  of  algebraic  equations  in  A,  B,  C: 


(61) 


(CiX^  +  C2)A  -  XB  -  XC4C  =  0, 


T^  +  (CrX^  +  c  )B  +  (cg  -  r^X^)C  «  0, 


C4rA  +  fc^  -  C2X^)B  (c^X^  Cjg)C  »  0. 


(62) 
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j;ii«Uiiiu.Aiu»i<'i'Rwy»KtMwm«awwa>a«a5Cwew»ii!»^^ 


This  system  of  equations  will  have  a  nontrivial  solution  provided  that 
the  determinant  of  the  coefficient  matrix  vanishes.  Expansion  of  the 
determinant  leads  to  the  following  polynomial  of  degree  six  in  T 
(bicubic  in  '^) : 

CiCCgCg  -  C^)^^  +  [CjCgCjQ  ♦  Cq  (c^Cg  +  CjC^)  +  2c^c^  + 

*  ‘^^5  -  *  '=2‘=6'=9  *  “=10  '‘=2^  *  ' 

*  2‘'4‘=S  *  “=^6  ■  =1=8  *  *  Cj„)P  *  C2<'=6'=10  '  '  "• 

C63) 

But  with  w=0,  (63)  reduces  to 

T^{cj(c5Cg  -  •*■  [cjCcgCjQ  ♦  CgCg  ♦  2c^Cg)  ♦  2c^Cy  * 

*  '^5  *  =91’'^  *  'd'6'io  -  'i>  •  *  '=^6  ’  'lo'  ■  “• 

C64) 

A  numerical  study  of  (64)  has  revealed  that,  for  the  range  of 
interest  here  (numerical  values  for  the  laminate’s  parameters  are 
given  in  the  next  section),  the  roots  of  (64)  may  be  expressed  as 

T  =  0,  0,  i  inj  i  n2> 

^2  >  0,  i  =  (-1)^^^ 
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Consequently »  the  functions  stated  in  (61)  now  assume  the  following 
forms : 


wCx) 

♦  (x) 


Aj  ♦  A^x  ♦ 
Bj  ♦  B2X  ♦ 
Cj  +  C2X  ♦ 


Ajcosn^x  ♦ 
Bjcosnjx 
Cjcosnjx  ♦ 


A.sinn.x  ♦ 
4  1 

B^sinOjX  + 
C^sinnjX  ♦ 


AjCoshn2X  ♦ 
B5coshn2X  ♦ 
Cj.coshn2X  ♦ 


Agsinhn2X, 

BgSinhr}2X, 

Cgsinhn2X. 

(65) 


The  eighteen  constants  A^,  Bj,  Cy  j  •  1»2,3,4,5,6,  in  (65)  are 
interrelated.  Substitution  of  (65)  into  (60)  yields  a  set  of  a  dozen 
relationships  among  those  coefficients.  It  is  a  straightforward 
algebraic  exercise  to  demonstrate  then  that 


S' 

R3C,. 

A3  «  RjC^, 

A4  »  -  RjCj, 

S  " 

S' 

RjCj, 

“  Vl* 

h  • 

Bj  ■  -  SjCj, 

h  ’ 

-  SjC^, 

85  ■  - 

C2  -  0, 

where 


*^1  " 

“  (OiMi)[V^lV^lVl^* 

A^»ni[l*c^(c,  -c^n^)]. 


^2^tc4C^-Cg+(c7^c^C5)n2]. 

(O7/A2)  ’ 

^2(1  ♦  Cj(Cg  ♦  CgO^)], 


(c 


10 
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Consequently,  (65)  may  now  be  expressed  as 


w(x)  =  A  +  R_C,x  +  R.C.cosn.x  -  R,C_sirtr),x  ♦  R^C^-coshn-jX  ♦ 
1  3  1  14  1  13  1  26  2 


♦  R^CgSinhn^x, 


il'(x)  =  SjCj  -  SjCjCosn^x  -  SjC^sinn^x  -  S2CgCoshn2X  -  S2CgSinhn2X, 
‘J'(x)  =  ♦  Cjcosn^x  ♦  C^sinHjX  +  CgC0shn2X  +  C^sinhn2X. 


If  both  ends  cf  the  beam  are  clamped,  then,  in  view  of  (49)  the 
boundary  conditions  are 


w(0)  =  <i-(0)  »  «^(0)  =  w(l)  *  1/^(1)  =  4'(1)  »  0. 


Inserting  (66)  into  the  first  three  boundary  conditions  in  (67),  we 


obtain 


Ai  +  RjC^  +  ^2^6  “ 


^3^1  ■  "  Vs  " 


Cj  ^  ^3  ^  ^5  * 


it'i  thf’  help  of  these  relationships,  we  can  write  (66)  as 


w(x)  =  Rj(cosnjX  -  1)C^  +  (*^3^21^  * 


R2sinhn2X)C5  +  R2(coshn2X  - 
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i|<(x)  «  -  SjC^sinnjX  ♦  (S3521  ♦  Sj532CO*niX  -  S2COshn2x)Cg  - 
-  S2CgSinhn2X« 


♦  (x)  ■  C^sinOjX  ♦  (^21  -  532CosnjX  coshn2X)C5  ♦  CgSinhn2X» 


with 


Si  '  (S  -  S2  *  IS3  *  SjVCs,  ♦  Sj), 

Substitution  of  (68)  into  the  last  three  boundary  conditions  in 
(67)  leads  to  the  following  homogeneous  system  of  algebraic  equations 
in  C^,  Cg,  and  Cgt 


^ij  ^j+3  * 


(69) 


where 


All*  R2(cosn2  -  1),  ■  ^3^21  *  *^1^32®^*'^! 

Ai3  »  R2(coShn2  -  1)»  A2J  *  -  SjSinnj* 

^2  “  ^3^21  ^1^32*^®*'^!  "  ’^2®®*^'^2*  ^23  "  “  S2Sinhn2» 

Ajj  »  sinnp  A22  ■  ^21  "  coshn2» 

Ajg  »  sinhn2» 


The  system  of  equations  (69)  will  have  a  nontrivial  solution  if  and 
only  if 

Det  (A.j)  ■  0, 

^  (70) 
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Expansion  of  the  detenninant  in  (70)  yields  the  relatively  simple 
relationship 

2R2(coshn2  -  l)sinnj  ♦  * 

♦  sinhn2  “ 

However,  if  we  introduce  the  identities 

sinrij  *  2sin(nj/2)cos(nj/2) ,  sinhn2  “  2sinh  (112/2) cosh  (1)2/2) , 

l“Cosr)j  *  23in^(nj/2),  coshn2  -  1  =  sinh^(n2/2), 

into  (71),  we  obtain,  after  some  rearrangement, 

sin(nj/2)  sinh(n2/2)  [2K^snh{r\2/2)  +  2Rj5j2^®*^Cnj/2)  + 

*¥21!  =■  o- 

As  in  the  case  of  the  Timoshenko  theory  (see  the  Appendix)  the 
smallest  positive  value  of  can  be  shown  to  arise  from 

sin  (n^/2)  »  0, 

—2  2  2 

Thus,  ®  2Tr,  and  since  in  (64) A  «  -  ~  (2Tf)  »  we  find  that  the 

value  of  A^j.  is  to  be  computed  from 

Cj(c5C9-c2)(2n)4  -  [^^(CgCj^  +  c^c^  +  2c^z^)  *  2c^c^  +  c^c^  + 

+  Cg](2n)2  +  c^(c^CjQ  -  cp  -  2c^Cg  ♦  c^c^  +  Cjg  «  0, 

or  the  expanded  form  of  this  result  which  can  be  shown  to  be  a  poly¬ 
nomial  of  degree  three  in  A,  A  discussion  of  the  numerical  results 
obtained  from  (72)  will  be  given  in  the  next  section. 
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Let  us  next  address  the  vibration  problem  for  a  clanqped-clamped 
laminated  beam.  For  we  must  use  (63)  to  compute  the  character¬ 
istic  exponents  T,  Considering  (63)  as  a  cubic  in  we  have  been 
able  to  show  that,  for  the  range  of  interest  here,  (63)  has  either 
one  negative  root  and  two  positive  roots  or  two  negative  roots  and 
one  positive  root.  For  the  present  we  confine  our  attention  to  the 
case 


T®  >  >  0. 

5  2 


Now  let 


"2  *  ■^2’ 


3  3’ 


it  being  understood  that  n^>0,  j  «  1,2,3.  Therefore,  (61)  leads  to 

w(x)  »  AjCoshjX  ♦  A2SinniX  ♦  AjCoshn2X  ♦  A4sinhn2X  + 

+  Agcoshn^x  A^sinhOjX, 

>l»(x)  =  BjCosnj,x  B2SinnjX  BjCoshn2  +  B^sinhn2X  ♦ 

+  BgCoshOjX  +  Bgsinhn3X, 

'{>(x)  »  Cjcosnix  +  C2sinnjX  +  C3Coshn2X  ♦  C^sinhn2X 
♦  C5Coshn2X  +  Cgsinhn3X. 


Proceeding  as  we  did  with  (65),  we  can  show  that 


Ai  a  R2C2, 

A2  »  -  RiCj, 

A3  " 

^4  *  ^2^ 

^5  "  *^3^6 » 

A^  “  R3C5, 

h  •  SjCj, 

B2  «  SjC, 

B4  =  s^c^. 

B5  »  S3C5, 

^6  “  S3C, 
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•'•-•Mvawi 


Kn»M^MiiwgBdaa»^  i3S*«aai^>i%«afff»^aa<^^ 


where 


\  “  Cni/Aj)  [c^c^  -  Cg  -  (C4C5  + 

R2  »  (^2/^2^  ‘  H* 

R3  =  (nj/Aj)  [c^cg  -  eg  •»•  IC4CS  ♦  c,)n3]. 

51  =  4“^  fcicynj  +  (C4  ♦  CjCg  -  C2C^)nj  -  c^Cg], 

52  =  [cjc^n^  -  (C4  ♦  CjCg  -  C2C7)n^  -  c^Cg], 

53  =  A3I  [cjc^n^  -  (C4  *  CjCg  -  c^cpn^  -  c^Cg], 


Aj  =  CjCgnJ  -  (1  +  CjCg  ♦  *  Ye* 

j  O 

S  “  Ys^2  *  *  ^1*^6  *  '*'  ‘^2*^6* 


A  O 


Therefore,  (73)  becomes 

w(x)  a  Rj[C2Cosn^x  -  CjSinHjX]  +  R2[c4Coshn2X  +  ( 
+  Rj(CgCoshn3X  +  C3Sinhn3x], 

ij^(x)  =  Sj[CjCosn^x  +  C2sinr’jX]  +  [C3Coshn2X  +  ( 
+  S3fC3C0stjn3X  +  CgSinhn3x], 

(|)(x)  =  CjCOSMjX  +  C2sinii2X  +  C3COShn2  +  C4Sinhn23 
•♦■  C3Coshn3X  +  Cgsinhn3X, 


;3sinhn2x]  + 


:4sinhn2x]  + 


(74) 
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Invoking  the  first  three  botmdary  conditions  in  (67),  we  can 
reduce  the  number  of  constants  in  (74)  to  three: 

w(x)  »  R2R2(coshri2X  -  cosnjx)C*  ♦  [R2(S2  -  S2)sinnix  + 

♦  R2(S3  -  Sj)sinhn2X  ♦  (Sj  -  S2)R3Sinhn3X]C5  + 

♦  R2R3(coshn3X  -  cosniX)Cg, 

i/'(x)  =  (S2RjSinhn2X  -  SjR2SinnjX)C^  +  [52(82  -  S3)cosn2X  + 

♦  32(8 j  -  S2)coshn2X  +  S3 (52  “  52)coshn3x]C5  + 

♦  (SjRjSinhrijX  -  SjR3sinn2x)Cg, 

<{i(x)  =  (RjSinhn2X  -  R2sinn2X)C^  +  [(52  -  52)cosn2X  + 

+  (S3  -  S2)coshn2X  ($2  -  S2)coshn3x]C*  + 

♦  (RjSinhn3X  -  R3sinn2X)C^, 


where  we  have  set 


R1C4, 


The  frequency  equation  is  now  derived  by  inserting  (75)  into  the 
last  three  boundary  conditions  in  (67),  In  the  usual  manner,  we 
again  obtain  (70),  where  now 

All  *  “  cosn2), 

Aj^2  “  R2^®3"^1^®^”^'^2  *  (S2"S2)R3Sinhn2, 

Ai3  *  RiR3[^oshn2  -  cosn2),  A22  *  S2RjSinhn2  -  SjR2sinnj, 

^22  **  52  (S2~S2)cosn2  ♦  S2(S2-Sj)coshn2  +  82(52-52) coshn3. 
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^23  *  ’^3  ’  S^R^sintij,  Ajj  »  RjSinhn2  -  R2sinnj, 

Aj2  *  (S2-Sj)cosni  ♦  (S2-Sj)coshn2  +  (5^-82) cosh  113, 

A33  =  RjSinhn3  "  RjSintij. 


Hxpansion  of  the  determinant  in  (70)  leads  to  the  frequency  equation 


2R1R2Y23Y  jjCi-cosHj  coshn2)sinhnj  +  2R2R3Yj2'^31^^"^°®^'^2 
coshn2)sinn2  ♦  2RjR3Y22Y23(l“Cosnj  coshn3)sinhn2  + 

+  (RjY^  -  *^3Yi2  "  sinhn2  sinhUj  =*  0, 


where 


^23  “  ^2  ’  ^3*  “  ^2* 


31  -  <^3  ■  '^1' 


^12  “  “"1  ’  '^2' 


2  2 

If  the  value  of  0  is  such  that  n2  assumes  a  negative  value, 
2  —2 

then  we  set  ^2  “  ”  *^2*  identities 


sinh  (iru)  =  i  sin  n  ,  cosh  (i~  )  -•  cos  tL, 


we  can  now  express  (76)  in  the  form 


2R1R2T23Y31  Cl”COSrii  cosTi2)sinhn3  +  2R2R3Y22731  (l“'^0S’i2 


coshn3)sinni  +  ^^1*^3~12~  '^l-cosrij  coshn3)sinn’2  + 
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where 


\  °  *  'l°6  *  *  °2'6’ 

^2  ■  ^2^^2>  f'4‘'-6  -  '8  -  *  '75^1  • 

^2  ■  ’^2"'  *  ^"4  ’  ‘I'S  •  ''2‘7j''b  -  ‘=2'=8J> 

^23  =  ^2  ‘  ^3>  ’'12  ■  '  ^2’ 

2  2 

For  a  certain  set;  of  values  of  u  ,  say  w  ,  n  =  1,?,3,.,,,  the 

n 

left  side  of  (76  or  77)  w;lll  vanish.  In  the  numerical  work  we  assign 
2 

a  value  to  u  and  compute  the  nj’s  (i.e.,  the  X^’s)  from  (63),  The 
2 

value  of  u  and  are  then  inserted  into  (76  or  77),  We  evaluate 

the  left  side  of  (76  or  77),  and  if  its  value  is  differt  jt  from  zero 

2 

a  new  value  is  assigned  to  u  and  the  process  is  repeated  until  a 

2 

value  of  u  that  causes  the  left  side  of  (76  or  77)  to  vanish  is 
found. 

NUMERICAL  RESULTS 

Sun  [1]  compared  the  dispersion  curves  obtained  fro.'  the  effective 
modulus  theory  and  the  microstructure  theory  for  laminated  beams,  using 
the  following  numerical  values  for  the  parameters; 


Y  “ 

100, 

0  »  2, 

n  =  0,8, 

5  = 

• 

CO 

Vj  »  0,2, 

=  0,35 
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Tne  values  of  n  and  C  given  in  (78)  correspond  to  a  composite  beam 
consisting  of  five  stiff  layers  and  four  soft  layers.  We  shall  also 
use  these  values  here  as  well  as  other  values  of  5,  i.e.,  the  number 
of  layer  pairs  in  the  composite  will  be  varied. 

The  numerical  results  obtained  here  from  the  microstructure 
theory  are  compared  with  the  corresponding  results  obtained  from 
the  effective  modulus  theory  (see  the  Appendix) ,  which  is  comprised 
of  Timoshenko  beam  theory,  with  effective  moduli  and  density, 
including  the  effect  of  a  uniform  axial  compressive  load.  For  a 
two-  phase  composite  medium  of  the  type  under  consideration  in  this 
investigation,  Sun  [1]  assumed  the  effective  Young’s  modulus, 
effective  shear  modulus,  and  effective  mass  density  to  be 

E  =  nE^  +  Ci'^n)E2,  u  =  nuj  +  (l“n)u.,, 

p  =  no,  +  (l-n)p,, 

(79) 

respectively. 

rhe  iiinged-Hinged  Beam 

Tlie  critical  value  of  the  buckling  parameter  X  and  the  natural 
fro:[uencies  of  vibration  are  computed  from  (57)  for  a  hinged-hinged 
beam.  For  a  compressive  load  P  applied  in  a  conservative  manner,  we 
may  calculate  the  value  of  X^^  from  (57)  by  setting  the  frequency 
parameter  u  equal  to  zero  and  n=l  in  (56). 
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In  Figure  2,  the  variation  of  the  critical  load  coefficient, 

X  ,  is  plotted  against  the  depth-to-length  ratio,  c»  over  the  range 
0<i;<0,2,  The  corresponding  results  obtained  from  the  effective 
modulus  theory  (see  the  Appendix)  are  also  shovm.  The  microstructure 
and  the  effective  modulus  theories  are  in  rather  close  agreement  for 
relatively  long  beams.  However,  as  the  values  of  c  increases,  the 
values  of  X^^.  as  calculated  from  the  effective  modulus  theory  decrease 
relatively  slowly,  whereas  the  values  as  computed  from  the  microstruc¬ 
ture  theory  decrease  very  rapidly.  At  C  »  C.2,  the  microstructure 
value  is  approximately  38%  of  the  effective  modulus  value.  The  effect 
of  increasing  the  number  of  layer  pairs  is  also  shown  in  Figure  2. 
Increasing  the  number  of  layer  pairs  beyond  twenty  has  a  negligible 
effect  on  the  value  of  X  .  The  values  of  X  for  the  beam  consisting 
of  only  five  stiff  layers  and  four  soft  layers  are  only  a  few  percent 
greater  than  those  for  a  beam  consisting  of  a  large  number  of  layer 
pairs. 

Eq,  (57)  may  be  considered  as  a  polynomial  of  degree  three  in 

2 

u  ,  Hence,  for  each  choice  of  the  integer  n,  (57)  will  yield  three 
branches  on  a  frequency  plot,  and  these  branches  may  be  labelled  the 
(i)  flexural,  (ii)  thickness-shear,  and  (iii)  microstructure  rotational 
branches.  Since  this  third  branch  is  associated  with  relatively  high 
frequencies,  we  shall  concentrate  our  attention  in  this  study  on  the 
more  important  flexural  and  thickness-shear  branches.  In  Figures  3-5, 
the  variation  of  with  t,  for  the  first  three  flexural  modes  (n=l,2,3) 
with  X  a  -  5,0,5,  is  plotted  for  a  beam  made  of  five  stiff  layers  and 
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EFFECTIVE  MODULUS  THEORY 


5 


igure  2.  Variation  of  the  critical  load  coefficient  A.  with 
depth-to-length  ratio  C  for  a  hinged-hinged  ELain. 


four  soft  layers.  For  small  values  of  c,  i.e,,  relatively  long  beams, 
the  effective  stifftiess  and  microstructure  results  are  in  excellent 
agreement,  however  as  the  value  of  c  increases  the  microstructure  fre¬ 
quencies  are  considerably  less  than  the  effective  modulus  frequencies, 
the  magnitude  of  the  effect  increasing  with  the  mode  number  n.  ‘s  is 
to  be  expected,  the  presence  of  a  tensile  load  (negative  X)  tends  t 
increase  the  natural  frequency,  whereas  a  compressive  load  (positive  A 
tends  to  decrease  it  relative  to  its  value  in  the  absence  of  an  axial 
load. 

In  Figures  6-8,  the  first  three  flexural  frequencies  for  X  «  0 
are  again  plotted,  however  now  the  effect  of  the  number  of  layer 
pairs  on  the  frequency  is  demonstrated.  An  increase  in  the  number 
of  layer  pairs  tends  to  decrease  the  value  of  the  difference 
becoming  more  pronounced  as  t  and  the  mode  number  n  increase.  For 
sufficiently  small  values  of  c  the  effect  of  increasing  the  number 
of  layer  pairs  is  virtually  negligible.  In  addition,  increasing 
the  number  of  layer  pairs  beyond  twenty  has  very  little  effect  upon 


the  values  of  the 


s 


Figure  3,  Variation  of  the  fifst  flexural  frequency  ui(n=l)  with  the 
applied  load  parameter  X  and  the  depth-to-length  ratio  ^ 
for  a  hinged-hinged  beam. 
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rigi  re  5.  Variation  of  the  third  flexura.l  frequency  w(n=3)  with  X 
and  C  for  a  hinged-hinged  beam. 


rigure  6.  Variation  of  the  first  flexural  frequency  w(n=l)  with  ^ 
and  the  number  of  layer  pairs  for  a  hinged-hinged  beam. 
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Variation  of  the  second  flexural  frequency  w(n=2)  with  ^ 
and  the  number  of  layer  p'lirs  for  a  hinged-hinged  beam. 
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Figure  8 


Variation  of  the  third  flexural  frequency  w(n=3)  with  i, 
and  the  number  of  layer  pairs  for  a  hinged-hinged  beam. 


In  Figure  9,  the  lowest  thickness-shear  frequency  is  plotted 
as  a  function  of  c.  In  the  range  |xl<5,  the  calculations  performed 
here  were  virtually  insensitive  to  the  values  of  X,  The  effective 
modulus  theory  leads  to  values  for  that  are  again  greater  than 
the  thickness-shear  frequencies  computed  from  the  microstructure 
theory.  Based  upon  the  microstructure  theory,  the  frequencies 
for  the  first  three  thickness-shear  modes  were  determined,  and 
their  variation  with  c  is  shown  in  Figure  10. 
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Figure  10.  Variation  of  the  first  three  thickness-shear  frequencies 
with  5  as  determined  from  the  microstructure  theory  for 
a  hinged-hinged  beam. 


The  Clamped-Claroped  Beam 

The  critical  values  of  the  buckling  coefficient  X  are  calculated 

from  (72).  In  Figure  11,  the  variation  of  X  is  plotted  as  a  func- 

cr 

tion  of  c  over  the  range  0<c<0,2,  Results  of  the  microstructure 
calculations,  (72),  and  the  effective  modulus  calculations  (see  the 
Appendix)  are  shovim.  For  very  long  beams  both  theories  predict 
essentially  the  same  values  for  X^^,  but  as  c  increases,  i.e.,  as 
shorter  and  shorter  beams  are  considered,  the  microstructure  theory 
predicts  much  lower  values  for  the  critical  load  coefficient  than 
does  the  effective  modulus  theory.  In  fact  at  c  *  0.2,  the  micro- 
structure  value  for  X^j.  is  approximately  one-sixth  of  the  effective 
modulus  value. 

Using  the  frequency  equations  (76)  and  (77)  in  conjunction  with 
the  characteristic  polynomial  (63),  we  have  computed  the  frequency 
parameters  a»j^,  n=»l,2,3,  for  the  lowest  three  modes  in  a  clamped- 
clamped  beam.  These  retults  are  shown  in  Figures  12,  13,  and  14, 

It  is  evident  that  the  values  obtained  from  the  effective  modulus 
and  microstructure  theories  agree  only  for  small  values  of  5,  For 
C>0,02  the  respective  curves  diverge,  and  the  values  of  predicted 
by  the  microstructure  theory  are  considerably  less  than  the  correspond¬ 
ing  values  computed  from  the  effective  moduli^  theory.  In  order  to 
show  the  effect  of  initial  tensile  or  compressive  axial  loads,  we  have 
plotted  the  values  of  as  factions  of  c  for  X  »  -  150,0,150  for  the 
effective  modulus  theory  and  for  X  =  -  50,0,50  for  the  mi crostructui 2 
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EFFECTIVE  MODULUS  THEORY- 


s 


Figure  H.  Variation  of  the  critical  load  coefficient  X  with  the 
depth-to-length  ratio  5  for  a  clamped-clampe8^beam. 


igure  12. 


Variation  of  with  ?  and  X  for  a  clamped  ■•clamped  beam. 


theory.  Indeed,  it  should  bo  obsor/ed  in  Figure  12  that  for  A  =  50 

the  value  of  the  frequency  parameter  Uj  starts  at  zero  for  ?  =  0  and 

increases  to  a  maximum  as  c  increases  and  then  decreases  to  zero  as 
i;  increases  to  C  =  0,1945,  Since  =  0  for  A  s  50  and  r,  =  0,1945, 
this  value  of  C  corresponds  to  the  value  ol  the  critical  buckling 
coefficient  for  a  clamped-c lamped  beam,  the  length-to-depth  ratio 
of  which  is  0,1945  (compare  Figure  11).  The  curves  for  and  uj 

with  A  =  so  show  similar  characteristics  but  approach  rero  at 

higher  values  of  c. 
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Figure  13.  Variation  of  with  i;  and  X  for  a  c.lamped-clamped  beam. 
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Figure  14.  Variation  of  u)^  with  C  and  X  for  a  clamped- clamped  beam. 


CONCLUSIONS 


Based  upon  the  numerical  results  reported  herein  for  the  buckling 
coefficients  and  natural  frequencies  of  vibration  of  a  laminated  beam 
subject  to  either  hinged-hinged  or  clamped-clamped  end  conditions,  it 
is  evident  that  the  effective  modulus  theory,  i.e,,  Timoshenko  beam 
theory  for  a  transversely  isotropic  bar  subject  to  initial  axial  stress, 
appears  to  be  reliable  for  engineering  design  purposes  only  in  the  limit 
of  fairly  long  beams.  The  extended  version  of  the  microstructure  theory 
that  was  developed  and  investigated  in  the  present  study  predicts  much 
lower  buckling  coefficients  and  natural  frequencies  in  laminated  beams 
with  depth-to- length  ratios  of  moderate  size,  the  magnitude  of  the 
effect  being  more  pronounced  in  the  case  of  the  clamped-clamped  beam. 

In  view  of  the  dispersion  curves  reported  by  Sun  [1],  this  was  to  be 
expected. 

Based  upon  the  effective  modulus  theory,  the  results  reported  in 
References  [6]- [9]  for  buckling  coefficients  and  natural  frequencies 
indicated  much  lower  values  than  those  predicted  by  the  classical 
Ruler- Bernoulli  beam  theory,  which  does  not  account  for  the  disparity 
in  the  longitudinal  and  transverse  elastic  moduli.  Indeed,  the  results 
become  more  and  more  divergent  as  the  ratio  of  the  transverse  Young’s 
modulus  to  the  longitudinal  shear  modulus,  D/G,  increases,  the  effects 
being  particularly  pronounced  when  E/G  exceeds,  circa,  50.  The  numeri¬ 
cal  results  presented  in  this  report  indicate  still  lower  values  for 
the  buckling  coefficients  and  the  natural  frequencies,  even  though  the 
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E/G  ratio  was  only  slightly  greater  than  3.  Consequently,  the  validity 
of  adopting  a  transversely  isotropic  beam  to  model  a  laminated  beam 
must  be  seriously  questioned  and  challenged. 

It  is  currently  being  proposed  at  the  Watervliet  Arsenal  that  an 
experimental  study  of  the  buckling  and  vibrational  response  of 
laminated  beams  and  plates  be  initiated.  Such  a  study  will  provide  a 
basis  for  rendering  further  judgment  on  the  ranges  of  applicability  of 
both  the  microstructure  and  effective  modulus  theories. 
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APPENDIX 


EFFECTIVE  MODULUS  THEORY 

With  the  aid  of  (16)- (19),  the  fimdamental  differential  equations 
for  the  effective  modulus  theory  may  be  derived  in  a  very  simple 
'ashion.  In  view  of  the  studies  reported  in  [6] -[9],  it  is  known 
that  the  QXtensional  deflection  u(x^,t)  in  the  beam  is  uncoupled 
from  the  transverse  deflection  w(Xj^,t)  and  the  rotation  (j>(Xj,t). 

Hence,  the  function  u(Xj,t)  may  be  neglected.  Thus,  we  replace 
(16)-(18)  by 

U  =  i  [El4.^j  ♦  kpA(w^j  -  (^)^], 

(A-1) 

T  =  +  Aw^], 

^  (A-2) 

n  «  'T  Oil  / 

^  (A-3) 

respectively. 

Therefore,  application  of  (19)  in  the  familiar  manner  leads  to 
the  following  equations  of  motion: 

<nbh  (w  -<}i,)-Pw,,-  pbhw  =  0 , 

(A-4) 

(E-P/A)I.(;i  ,,  +  kijA(w  .  -  (j))  -  pi.  4.  =  0, 

’  '  (A-5) 
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with  the  boundary  conditions 


(;.)  either  tewACw  ,-^)+Pw,=0  or  6w  =*  0 

•  *  ♦  • 

(ii)  either  (E-P/A)l^(l)  ^  *  0  or  6ij(  ■  0 

CA-6) 

on  Xj  «  0,*..  We  have  set  »  P/A, 

—  2  i/2 

Setting  Xj  »  lx,  (j)  =  ({i/l,  and  t  =  (P2*-  /U2)  "f*  oxpress 

CA-4)  and  (A-5)  as  follows: 

(b,  -  Xb2)w”(x,T)  -  bi?*(x»T)  »  b,.w(x,T), 

^  (A-7) 

(btn  -  Xb*  )^'*(x,t)  ♦  b,(w'  -  «  b*  ?(x,r), 

^  1®  (A-8) 

where  bj,  b2,  and  bj3  are  given  in  the  text  and 

^  6c2(l+ed)  ^  ^  ^  1+ed 

12k  (1+d)  144k  12k  1+d 

For  simply-supported  ends,  the  boundary  conditions  for  the  beam 
are  w  =  p  =  0  at  x  =  0,1,  If  we  insert 

w(x,t)  =  A^sin  mrx  cos  ojt,  TCXjt)  »  A^cos  nnx  cos  cot 

into  CA-7)  and  (A-8),  we  obtain,  after  some  manipulation,  the 
frequency  equation 

-  t^bi3.(nn)2  fb^^bj^*  b^bj^-  ^(bisb^i^b^bj^)]  . 

+  (nTi)^[(n’T)^CbjQ-  AbJpCbj-  Xb2)  -  Xbjb2)  »  0 

(A-9) 
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The  frequencies  for  the  flexural  and  thickness-shear  modes  may  be 
obtained  from  (A-9)  with  the  aid  of  the  quadratic  formula.  Numerical 
results  calculated  in  this  manner  have  been  plotted  in  Figures  3-5 
and  9. 

In  the  case  of  buckling,  the  value  of  the  critical  load 
coefficient  is  obtained  by  setting  w-O  and  n«l  in  (A-9) : 


^^(bJo-  ^bJpCbi-  Xbp  -  Xbib2  -  0. 


It  is  easily  shown  that 


2”^l>nl>2  *cr  ■  - 


(A-10) 

Numerical  results  derived  from  (A-10)  have  been  plotted  in 
Figure  2. 

Let  us  consider  the  buckling  and  vibration  of  a  clamped-clamped 
beam  separately.  In  the  former  we  may  assume  that  w  and  ^  are 
independent  of  time  t,  so  that  (A-7)  and  (A-8)  reduce  to  the  follow¬ 
ing  pair  of  ordinary  differential  equations: 


(bj-  Xb2)w‘'(x)  -  bj?'(x)  3  0, 


(A-11) 


it 


(bjo”  Xbn).fr"(x)  ♦  bj  [W(x)  -  Hx)]  -  0. 


(A-12) 
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The  solutions  of  (A-11)  and  (A-12)  are  easily  found  to  be 

TCx)  =  A,  cosAx  ♦  A  sinAx  +  A3, 

^  (A-13) 

w(x)  =  A4  +  A-x  +  Cc/A)A,  sinAx  -  (c/A)A~  ccsAx, 

(A- 14) 

where 

c  =  bj/(bj  -  Ab2),  A^  *  Abjb2/(bj  -  Ab^)  (b*^-  AbJj). 

The  boundary  conditions  for  a  clampod-clamped  beam  are 
T=w  =  0  at  X*  0,1. 

(A- 15) 

Substitution  of  (A-13)  and  (A-14)  into  (A-25)  leads  to  a  system  of 
homogeneous  algebraic  equations  in  the  i  »  1,2, 3, 4,  which,  in 
the  usual  fashion,  leads  to  the  equation 

sinA  =  2(c/A)  (1-cosA) , 

(A- lb) 

from  which  we  can  calculate  the  value  of  the  critical  load  coefficient. 
Using  the  trigonometric  identities 

sinA  ®  2  sin  {A/2)  cos  (A/2),  cosA  «  1  -  2  sin^(A/2), 
ws  can  verify  that  (A-16)  leads  to 

sin  (A/2)  =■  0  and  tan  CA/2)  =  A/2c, 
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The  smallest  value  for  is  found  to  originate  from  sin  (A/2)  *  0, 

i,e,,  A  =  2ir,  Now  in  view  of  the  definition  of  A  above,  it  is  a 

straightforward  matter  to  demonstrate  that  A  is  to  be  computed 

cr 


from 


bzb*!  A^^  -  (bjbJi  +  b2bjQ  +  bjb2/4ir2)A_  +  b,b* 


cr  ‘"rio 


0, 


Numerical  results  obtained  from  this  expression  have  been  plotted 
in  Figure  11, 

To  determine  the  frequencies  of  vibration  of  a  clamped-clamped 
beam,  we  consider  (A-7)  and  (A-8)  subject  to  (A-15),  Inserting 

wCx,T)  *  WCX)  cos  OJT,  T(X,t)  =  (j>Cx)  cos  WT 

into  (A-7),  (A-8),  we  have 


(bj  -  Ab2)w*’(x)  -  bji{i’(x)  =  -  bj2W^w(x), 


(A-17) 


(bJo  -  ^bjj)4i”(x)  +  bj  [w»(x)  -  '))(x))  =  -  bjgU)^4(x), 


(A-18) 


If  we  introduce  a  new  dependent  variable  ij'(x),  according  to 
"{•(x)  =»  '//'(x), 

then,  without  loss  of  generality,  (A-18)  can  be  replaced  by 
(bjo  -  (x)  ♦  b^  [w(x)  -  iJ<(x)]  »  ..  bjgW^Kx), 


whence 


*  2, 


b^w(x)  -  (bj  -  bjgw  )ij»(x)  -  Cojo  -  *  (x). 


(A- 19) 
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Substitution  of  (A-19)  into  (A-17)  leads  to 

-  p  i^(x)  -  0, 


CA-20) 


where 


P2  *  P2/P4  and  p  a  Po/p  . 


Pq  *  I*  P2  =  «2+Xb  b 

16  Djg 


P4  =  Cbi  -  AbpcbJo  -  XbJi),  n2  »  b;//bj. 


.2  .  V*  .2, 


However,  it  can  be  shown  that  (A-.20)  may  be  represented 

*  6l)Hx)  »  0,  D  d/dx, 

where,  with  a  =  1,2, 


(A-21) 


Ci^  +  C2(l-n  )  ..  (-1)  [c^'*  +  + 

7  27  i/2 

+  cj  (i-n^)^]  , 


(A-22) 


Cl  =  ('>iV»'>n)  '  =2  »  ibjb. 


bjb 

C-  s 

'  b*  * 
16 


^4  =  -  ^^3)  (b^^  .  Abjp, 


Cl.  =  Xb^2  *  (bJo  -  bj  -  Abj). 
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It  is  now  important  to  observe  that 


0  for  >  0 


and  thu: 


(A-23) 


'11)6  significance  of  the  sign  of  6^  in  the  two  frequency  ranges  rests 
in  the  form  of  the  solution  of  (A-21).  By  the  principle  of  super¬ 
position  in  the  theory  of  differential  equations,  we  can  replace 
(A-21)  by 


with 


(D^  +  62)^^  (X)  =  0,  (d2  +  62)'J'2W  “  0* 


v(x)  = 


(A-24) 


(A-25) 


It  is  evident  from  (A-23)  and  (A-24)  that  'l'j(x)  will  always  be  given 
by 

ij;j(x)  =  Aj  COS  o^x  Bj  Sill  6jX, 


whereas 

ij'‘^(x)  =  A2  cosh  T2X  +  sinh 

with  ^  "  ^2’  or  the  form 

^^(x)  =  A2  cos  +  B2  sin  (S^x, 


(A-26) 


0  <  n  <  1, 


1  <  n. 


(A-27) 


(A-28) 
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In  terms  of  our  present  notation,  the  boundary  conditions 
(A- 15)  can  now  be  expressed  as 


w  =  sO  at  x  =  0,l» 


(A-29) 


The  second  condition  in  (A-29),  in  view  of  fA-25),  leads  to 


I  ♦ 

’’  0  at  X  *»  0,1. 

X  M 


(A- 30) 


By  virtue  of  (A-24)  and  (A-25),  (A-19)  becomes 


W(x)  = 


.  J.  (bj 


Xbjj)  iJ^jCx)  + 


-  <|»^(x), 


(A-31) 


Thus,  the  first  boundary  conditions  in  (A-29)  become 


Fi^i  +  F2X>2  =  0 


X  =  0,1, 


(A-32) 


where 


F„=  1-1)2*  (s2/b,)(b;„  -  XbJj), 


Consequently,  we  have  established  that  the  boundary  conditions  for 


the  clamped-clampod  beam  may  be  expressed  as  (A-30)  and  (A-32), 
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Case  I, 


Then 


2  —0 
0  <  <  1.  In  this  case  *^2 


=  1  -  i-52  .  (62/bp  (b*^  -  Abjp. 


Substitution  of  (A-26)  and  (A-27)  into  (A-30)  and  (A-32)  leads  to  a 
system  of  homogeneous  algebraic  equations  in  the  coefficients  A^,  B^j^, 
ct  =  1,2,  from  which  we  can  derive  the  frequency  equation 

2  F  F  6  ”3”  (1-cos  5,  cosh  F  )  +  (F^  "5?  -  F?6^)sin  6,  sinh  T  =  0, 
1212  1  2  1221^  1  2 

(A-33) 

Case  II.  1  <  fl.  In  this  situation,  (A-27)  is  replaced  by 
(A-28),  and  the  analysis  of  the  previous  paragraph  is  repeated. 

Because  ^2  ~  "  *^2  ^'^2^^^*  ^2  ~  ^^^2'  ^ 

s  =  +  1.  Now 

F2  =  1  -  4  (62/bp  (b*Q  -  Xb*p 


and 


cosh  6^  =  cos  5, 

2  i 


sinh  6^  =  is  sin  6  , 
2  2 


It  is  next  easy  to  verify  that  (A-33)  becomes 

2F,F^6,6_(l-cos6,  cos6,)  -  ifhl  *  F262)sin6,  sin6^  =  0, 

i  c  i  ^  i  ^  i  Z  CL  1  2 

(A-34) 

for  1  <  n.  Upon  computing  the  value  of  fi  from  (A-33)  or  (A-34),  we 
can  determine  10  from 


*  .1/2 


(u  =  n 
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